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ABSTRACT 

We have used 4752 days of data collected by the Birmingham Solar- 
Oscillations Network (BiSON) to determine very precise oscillation frequencies of 
acoustic low-degree modes that probe the solar core. We compare the fine (small 
frequency) spacings and frequency separation ratios formed from these data with 
those of different solar models. We find that models constructed with low metal- 
licity are incompatible with the observations. The results provide strong support 
for lowering the theoretical uncertainties on the neutrino fluxes. These uncertain- 
ties had recently been raised due to the controversy over the solar abundances. 

Subject headings: Sun: helioseismology - Sun: interior - Sun: abundances - 
neutrinos 
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1. Introduction 

One of the key inputs used in constructing models of stars, including the Sun, is the 
heavy element abundance, Z, or the ratio of the abundance of heavy elements to that of 
hydrogen ZjX. Helioseismic studies of solar models constructed with the Grevesse & Sauval 
(1998; henceforth GS) abundance, and the even earher Grevesse & Noels (1993) abundances, 
show a remarkable agreement with the Sun (e.g., Ghristensen-Dalsgaard et al. 1996; Bahcall 
et al. 1997; Morel et al. 1999; Basu, Pinsonneault & Bahcall 2000, etc.). The situation has 
changed in the recent years with the claim by Asplund et al (2001), AUende Prieto, Lambert 
& Asplund (2001, 2002), Asplund (2004), Asplund et al. (2004), Asplund et al. (2005), that 
the oxygen, carbon and nitrogen abundance in the Sun needs to be reduced. The reduction 
in the oxugen abundance implies the reduction of the abundances of other elements such 
as Ne, Ar also. This reduces the solar heavy element abundance from the GS value of 
ZjX = 0.0229 to ZjX = 0.0176. Asplund, Grevesse & Sauval (2005; henceforth AGS) 
presented a new table of abundances based on the above reduction. In addition, based on 
Asplund (2000) they reduced the sihcon abundance (and thus lowered the abundance of all 
the meteoritic elements), and this reduced ZjX further to 0.0166. Solar models constructed 
with the new, low abundances do not satisfy the helioseismic constraints of the position of 
the convection-zone base and the helium abundance in the convection zone. Models with the 
lower heavy-element abundances are found to have too shallow a convection zone (henceforth 
CZ), and a helium abundance, Y . that is too low. Furthermore, the sound-speed and density 
profiles of the models do not match those of the Sun either (Bahcall & Pinsonneault 2004; 
Basu & Antia 2004; Bahcall et al. 2005a, 2005b, 2005c, 2006; Delahaye & Pinsonneault 
2006) . By looking at the structures of the main heavy-element ionization zones in the solar 
convection zone Antia & Basu (2006) argue that the solar heavy element abundance has to 
be high, i.e., close to the old GS values. Pinsonneault & Delahaye (2006) by considering the 
physics of stellar atmospheres also conclude that high CNO abundances are favored. 

Almost all the helioseismic tests of solar models constructed with the new abundances 
have focused on the solar convection zone. Few studies have looked at what happens to the 
solar core. Where the core has been considered it has been through results of inversions, 
where inversion data on the core have been combined with results pertaining to the other 
parts of the Sun to give RMS deviations of the sound-speed and density profiles of the 
models with respect to those of the Sun (e.g., Bahcall et al. 2005a, 2006). In this work, we 
concentrate on the solar core. The solar core is, of course, where all the nuclear reactions 
take place, and where all the energy and neutrinos are released - hence the interest in what 
happens to the core when solar metallicity is reduced. 

Only modes of low degree (/ = 0-3) can be used to determine the structure of the solar 
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core. In this paper we report a study in which observed low-degree frequencies have been 
compared with those of several standard solar models constructed with different abundances. 
The solar p-mode data come from observations made by the ground-based Birmingham 
Solar-Oscillations Network (BiSON; Chaplin et al. 1996). The BiSON instruments make 
disc-averaged observations of the Sun in Doppler velocity. These observations have provided 
the field of helioseismology with high-quality, long-term monitoring of the low-/ modes. Here, 
we use frequencies determined to excellent precision from a BiSON timeseries lasting 13 years. 

The aim of this paper is to compare and contrast the frequencies of different solar models 
with those of the Sun in order to determine whether we can distinguish between models made 
with the older, higher abundances on the one hand and those made with the newer (AGS) 
abundances on the other. Our aim is then to see if we can rule out one set of abundances, 
thereby reducing the uncertainty in the solar abundances. This is particularly relevant for 
the study of solar neutrinos. Recently the theoretical uncertainty on solar neutrinos had 
been raised because of the uncertainty in the solar abundances (Bahcall & Serenelli 2005, 
Bahcall et al. 2006). The uncertainties in the values of the heavy-element abundances of 
the Sun are the largest source of the theoretical uncertainty in calculating the p-p, pep, ^B, 
13 ^^O, and ^''F solar neutrino fluxes (Bahcall & Serenelli 2005). The difference between 
the GS and AGS abundances is larger than the error-bars on each of the abundances, and 
hence Bahcall & Serenelli (2005) and Bahcall et al. (2006) calculated the uncertainties in the 
neutrino fluxes by designating the difference between GS and AGS values as the uncertainty 
in the composition. In this work we try to examine if data on the low-degree modes are of 
sufficient quality to distinguish between models with the older, higher abundances and those 
with the newer, lower abundances and hence to thereby reduce the neutrino uncertainties to 
pre- AGS levels. 

A simple comparison of the frequencies of the Sun and solar models is often not very 
informative. One reason is that all the modes sample the outer layers of the Sun, and hence 
even low-degree modes do not have localized information on the solar core. Another reason 
is that the physics used in the solar models does not reproduce the structure of the near- 
surface layers very well: the main culprit is the very basic treatment of convection, generally 
using the mixing length theory, which breaks down close to the surface, another factor is that 
the assumption of adiabaticity, used in the calculation of the frequencies of solar models, 
breaks down near the solar surface. Both these errors, and other uncertainties relating to the 
solar surface, give rise to frequency errors that are functions of frequency alone (e.g.. Cox & 
Kidman 1984; Balmforth 1992 and discussion in Christensen-Dalsgaard 2002). To study the 
core one therefore, normally compares the so-called fine or small frequency spacings of the 
low-Z p modes given by the combination. 

dii+2{n) = - i/„_i,i+2, (1) 
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where v is the frequency of a mode of degree I and radial order n. The frequencies Vn,i 
and z/„_i j_|_2 are very similar and hence are affected in a similar way by near-surface ef- 
fects. By taking this difference in frequency a large part of the effects from the near-surface 
uncertainties thus cancels out. 

The fine spacings are determined predominantly by the sound-speed gradient in the 
core. Using the asymptotic theory of p-modes it can be shown that (see e.g., Christensen- 
Dalsgaard & Berthomieu 1991) 

^lin) dcdr 

where, R is the solar radius, and A;(n) is the large frequency spacing given by 

A/(n) = - yn-i,h (3) 

which depends inversely on the sound-travel time between the centre and the surface of the 
Sun. Given that the gradient of the sound-speed is large in the core, and r there is small, 
the integral in Eq. 2 is dominated by conditions in the core, and hence the fine spacings are 
a useful tool to study conditions in the solar core. 

Although the difference between the frequencies of two modes that have very similar fre- 
quencies does reduce the effect of near-surface uncertainties, there are some residual effects. 
One way of reducing the effects of the near-surface errors is to use the so-called frequency sep- 
aration ratios. The frequency separation ratios (Roxburgh & Vorontsov 2003; Oti Floranes, 
Christensen-Dalsgaard & Thompson 2005; Roxburgh 2005) are formed from the fine (small) 
frequency spacings and large frequency spacings of the modes. The ratios are constructed 
according to: 

/ N do2{n) di3(n) 

Since both the small and large spacings are affected in the same manner by near-surface 
effects, these ratios are somewhat independent of the structure of the surface. 

In this paper we use data from BiSON to form observational estimates of the fine 
spacings, dQ'2.{n) and di3(n), and the separation ratios, rQ-2,{n) and ri3(n). The observed 
spacings and ratios have then been compared with estimates formed from a number of 
different solar models. BiSON data on the fine spacings have been used in several previous 
studies (e.g., Elsworth et al. 1990; Chaphn et al. 1997), but those here come from much 
longer observations and have far superior precision. 

The rest of the paper is organized as follows: the observed data are described in § 2, 
the solar models are described in § 2.2, we present and explain our results in § 3, and we 
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discuss the implications of these results to solar neutrino predictions in § 4, and we present 
our conclusions in § 5. 

2. Data 
2.1. Observations 

We have used Doppler velocity observations made by the BiSON over the 4752-d period 
beginning 1992 December 31, and ending 2006 January 3. Mode peaks in the power spectrum 
of the complete time series were fitted in the usual manner (e.g., see Chaplin et al. 1999) to 
yield estimates of the low-/ frequencies. 

2.1.1. Impact of the solar cycle on the Sun- as- a- star frequencies 

The BiSON observations are examples of the so-called Sun-as-a-star data: averages 
over the visible disc of the perturbations - here, in Doppler velocity - associated with the 
modes. Chaplin et al. (2005) showed that fine spacings and separation ratios made from 
the Sun-as-a-star data are sensitive to the changing surface activity along the solar cycle. In 
short, the acoustic asphericity leaves its imprint on the azimuthally dependent Sun-as-a-star 
frequencies. 

Some of the mode components are effectively missing from the data because the sensi- 
tivity of the observations to them is so low. This is a consequence of the visibility of any 
given m being a strong function of the angle of inclination, i, offered by the star. Extant 
Sun-as-a-star observations, such as those of the BiSON, are made from i 90 degrees. This 
means that only components with I -\- m even have non-negligible visibility. As such, it is 
not possible to estimate directly the frequency centroid; and estimates of the 'frequency' 
of a mode will be influenced strongly by the \m\ = I modes, which are most prominent in 
the data. The Sun-as-a-star flne spacings are therefore formed from modes with different 
combinations of / and m, which may suffer different sized frequency shifts through the solar 
cycle. 

For the do2{'n) data, estimates of the / = 2 frequencies are dominated by the \m\ = 2 
modes. They show significantly larger cycle shifts than the lone components of the nearby 
/ = modes (e.g., see Chaplin et al. 2004a; Jimenez- Reyes et al. 2004). The dQ2{n) therefore 
decrease in size as levels of surface activity increase. Differences between the |m| = 3 and 1 
components of, respectively, the I — 3 and 1 modes, are less pronounced, and so cycle shifts 
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arising in the dis{n) are smaller, at fixed frequency, than in their do2{n) counterparts. 

The separation ratios show a similar behavior. For these data, changes in the large 
spacing, A;(n), which lies in the denominator of Equation 4, are negligible. This is because 
A;(n) is formed from the difference of two frequency estimates in modes with the same 
combination of I and m. While it is true that there is a dependence of the mode shifts on 
frequency, the separation of the overtones used to form the A;(n) is sufficiently small that 
the impact of this dependence is modest. Changes to the separation ratios are therefore 
dominated by changes to the fine spacings. This means that fractional changes in the rQ2{n) 
are, to first order, the same as those in the do2{'n) [and likewise for the ri3(n) and di3{n)]. 



2.1.2. Correction of the Sun- as- a- star frequencies for the solar cycle 

Prior to calculation of the fine spacings and separation ratios, we therefore removed the 
solar-cycle shifts from the raw fitted low-/ frequencies. Our procedure rests on the assumption 
that variations in certain global solar activity indices can be used as a proxy for the low-/ 
frequency shifts, duneit). We assume the correction can be parameterized as a linear function 
of the chosen activity measure, A{t). When the 10.7-cm radio flux (Tapping & De Tracey 
1990) is chosen as the proxy, this assumption is found to be robust (e.g., Chaplin et al. 
2004a) at the level of precision of the data. 

Consider then the set of measured eigenfrequencies, I'neit), that we wish to correct, 
extracted from data collected over the t — 4752 d epoch when the mean level of the 10.7-cm 
radio flux was {A{t)) = 121 x 10~^^ Wm~^ Hz~^. We make the correction to the canonical 
quiet-Sun level of the radio flux, Aq^iet, which, from historical observations of the index, is 
usually flxed at 64 x 10~^^ Wm~^ Hz"^ (see Tapping & DeTracey 1990). The magnitude of 
the solar-cycle correction - which must be subtracted from the raw frequencies - will then 
be: 

Sl^neit) = 9e ■ J={V\ ■ {{A{t)) - ^quiet]- (5) 

The gi are /-dependent factors that calibrate the size of the shift: recall from Section 2.1.1 
above that the Sun-as-a-star shifts alter significantly with /. Owing to the nature of the 
Sun-as-a-star data, the change with / comes from the spatial dependence of the surface 
activity. To determine the (7/, we divided the 4752-day timeseries into 44 independent 108-day 
segments. The resulting ensemble was then analyzed, in the manner described by Chaplin 
et al. (2004a), to uncover the dependence of the solar-cycle frequency shifts on the 10.7-cm 
radio flux. The T\i'\ in Equation 5 is a function that allows for the dependence of the shift 
on mode frequency. Here, we used the determination of T\i'\ to be found in Chaplin et al. 
(2004a, b). 
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Uncertainty in the correction is dominated by the errors on the gi. These errors must 
be propagated, together with the formal uncertainties from the mode fitting procedure, to 
give uncertainties on the corrected frequencies, Vniii) — 5un£{t)- The corrected uncertainties 
are, on average, about 10 % larger than those in the raw, fitted frequencies. 

After application of the correction procedure, fine spacings and separation ratios were 
calculated according to Equations 1 and 4 respectively. Uncertainties on the corrected fre- 
quencies were propagated accordingly to give those in the spacings and ratios. 

Fig. 1 shows the differences between fine spacings made with and without the correction 
(in the sense corrected spacings minus raw spacings). The ratio data give similar-looking 
plots. As indicated previously, it is the dQ2{n) (and by implication the rQ2{'n)) that are most 
affected by the correction procedure. 

2.2. Models 

For the major part of the work we used eight solar models that have been constructed 
with different physical inputs. The models are: 

JCDS — Model S of Christensen-Dalsgaard et al. (1996). This model was constructed with 
surface Z/X = 0.0245 (Grevesse & Noels 1993), OPAL(1992) opacities (Rogers & Iglesias 
1992) and the OPAL(1996) equation of state (Rogers, Swenson & Iglesias 1996). This model 
was used because many helioseismological results in literature are based on this reference 
model. 

SAC — Model Seismici of Couvidat et al. (2003). It used OPAL(1996) opacity tables 
(Iglesias and Rogers 1996), and the OPAL(1996) equation of state (Rogers, Swenson & 
Iglesias 1996). The model has {Z/X) = 0.02628. 

BP04 - Model BP04(Garching) of Bahcall et al. (2005c). This model was constructed with 
GS abundances {Z/X = 0.0229), OPAL(1995) opacities (Iglesias & Rogers 1996) and the 
OPAL(2001) equation of state (Rogers 2001, Rogers & Nayfonov 2002). 

BP04H Model BP04+ of Bahcall et al. (2005a). This model was constructed with 

abundances from Asplund et al. (2000, 2004), Asplund (2000), AUende Prieto et al. ( 2001, 
2002). These abundances imply Z/X = 0.0176 as opposed to the Z/X = 0.0229 of GS. Only 
abundances of C, N, O, Ne and Ar are lowered in comparison to model BP04. The rest of 
the input physics is the same as BP04. 

BS05(OP) — Model BS05(OP) of Bahcall et al. (2005c). This model is similar to model 
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BP04, but has been constructed with opacities from the OP project (Badnell et al. 2005) 

instead of the OPAL opacities. There are some other subtle differences too, such as in the 
treatment of diffusion where, unhke in BP04, each metal has a different diffusion velocity as 
derived from Thoul et al. (1993). Also, unlike in BP04, ^^O is not burned at all. The model 
was constructed with GS abundances and has a surface Z/X = 0.02292. 

BS05(AGS, OP) — Model BS05(AGS, OP) of Bahcall ct al. (2005c). This model is similar 
to model BS05(OP), and was constructed with AGS abundances, having surface Z/X — 
0.01655. 

BBS05-3 — Model 3 of Bahcall et al. (2005b). This model has the same physics as BS05(OP) 
and BS05(AGS, OP), but the Ne, Ar, CNO and Si abundances have been enhanced with 
respect to the AGS values. The model has Z/X = 0.02069. 

S06+(AGS, OP) — A model constructed specifically for this work. It is similar to model 
BS05(AGS, OP), however, the opacities have been artificially increased by 13% in the tem- 
perature range of 2 to 5 million Kelvin to get the helioseismically determined position of the 
convection-zone base. 

We also used three other models to test the effect of heavy elements, keeping all other 
input physics the same. These models were constructed using YREC, the Yale Rotating 
Evolution Code, in its non-rotating configuration (Gucnther et al. 1992). All models have the 
same physics inputs except the heavy-clement abundances. The models use the OPAL2001 
equation of state, and OPAL (1996) opacities. The models are: 

YREG(AGS) — model with surface Z/X = 0.0165, the AGS value of surface metallicity. 

YREC(GS) — model with surface Z/X = 0.0229, the GS value of surface metaUicity. 

YREC(0.03) — model with a high metallicity, Z/X = 0.03, constructed to check if the effect 
of metallicity is monotonic. 

All models except S06+(AGS, OP) and the three YREC models arc published models. 
The models have been calibrated to slightly different values of luminosity and radius. The 
SAC model has a radius of 6.95936 x 10^° cm. All others have a radius of 6.9598 x 10^° cm. 
Model JCDS has a luminosity of 3.8456 x lO^^^ ergs/s; models BP04, BP04+, BS05(OP), 
BS05(AGS, OP), BBS05-3, and S06+(AGS, OP) have luminosities of 3.8418 x 10=^^ ergs/s; 
and the three YREC models have a luminosity of 3.851 x 10^^ ergs/s. The luminosity that 
the SAC model was calibrated to has not been published. 

The frequencies for each of these models was calculated by solving the full set of equa- 
tions describing stellar oscillations. The numerical precision of the frequencies was increased 
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using the Richardson extrapolation technique. Once the frequencies were calculated, The 
large and small separations, and the separation ratios were calculated using equations 1, 3, 
and 4. 



3. Results and Discussion 

Each panel of Fig. 2 shows the fine spacings, do2{ri), from one of the solar models. Model 
data are rendered as a sohd line, and the models are identified in the plot titles. The panels 

also show for direct comparison the dQ2{n) calculated from the corrected BiSON frequencies 
(points with error bars). Figs. 3, 4 and 5 show similar plots for the dis{n), ro2{n) and ri3(n) 
data respectively. 

The fine spacings and separation ratios show a marked dependence on mode frequency. 
Removal of these trends allows for a more detailed visual comparison of differences between 
the BiSON and model data. In the case of the fine spacings, we have followed an approach 
similar to that given in Chaplin et al. (1997). Here, a simple linear model of the form 

dii+2{n) ^ Co + ciu (6) 

was fitted to each set of BiSON fine spacings. The Fit to the do2{n) data, made over the 
range ~ 1408 to ~ 3985 //Hz, gave best-fitting coefficients of Cq = 15.95 ± 0.06 fiHz and 
Ci = (—2.24 ± 0.02) X 10"^; while the fit to the di3{n) data, made over the range ~ 1473 to 
~ 3640 /xHz, gave coefficients cq = 26.29 ± 0.20 /xHz and ci = (-3.35 ± 0.08) x 10"^. The 
best-fitting models were then subtracted from the BiSON spacings, and each set of model 
spacings, to yield the residuals that are plotted in the two panels of Figs. 6 (see caption 
for details). By adopting this approach of characterizing the spacings by a simple linear 
model we preserve in the residuals various features in frequency that are common to both 
the observations and models. The feature around ~ 2000 /xHz which is present in the various 
curves in Fig. 6 is the signature of the Hell ionization zone. The presence of this feature in 
the data shows that the fine-spacings are not immune to the near-surface structure. 

The separation ratios ro2(?^) in contrast show marked departures from linear behavior 
with frequency at the low-overtone end of the plotted data. Furthermore, these data, and 
the ri3(n), tend to vary more smoothly in frequency than do their fine-spacing counterparts. 
In Fig. 7 we therefore plot just the differences between the observed and model ratios (see 
caption for details). 

Tables 1 and 2 give quantitative measures of the differences between the observed and 
model datasets. Table 1 shows weighted mean differences (in the sense BiSON minus model). 



-10- 



and weighted RMS differences, for the fine spacing data (all in /xHz). The formal uncertain- 
tics on the BiSON spacings were used to fix the weights (with the usual uncertainty-squared 
Gaussian weighting applied), and to calculate an internal error on each mean or RMS differ- 
ence. Significance levels for the differences (in units of sigma) appear in the table in brackets, 
and were computed by dividing each mean or RMS difference by its associated internal error. 

Table 2 shows similar data for the separation ratios. Here, the mean differences are 
weighted fraciional differences between the BiSON and model values. Because each ratio is 
formed from the quotient of two separations in frequency, with the uncertainties on these 
separations assumed to follow Gaussian distributions, determination of the fractional dif- 
ferences allows errors to be given that are also Gaussian (and not asymmetric) in the final 
measure. 

From inspection of the various figures, particularly Figs. 6 and 7, and the data in the 
two tables, it is apparent immediately that the models that show the poorest correspon- 
dence with the observations are those that have lowered abundances, i.e., models BP04-I-, 
BP05(AGS,OP) and S06+(AGS,OP). Model BP04+ has shghtly higher metaUicity than 
models BP05(AGS,OP) and S06+(AGS,OP), and so fares shghtly better in the compar- 
isons. The models with higher abundances do much better: model JCDS, although con- 
structed with outdated inputs, agrees quite well with the observations, as do models BP04 
and BS05(OP). The non-standard model, BBS05-3, does well too. Model SAC docs not fare 
as well as models JCDS or BP04, however, it still agrees better than do the models with 
lower abundances. Thus we can say the models with lower than GS abundances have core 
structures that do not match that of the Sun. 

We shall not discuss models JCDS and SAC further, except to note again that they 
fare much better than the \ow-Z/ X models. Models JCDS and SAC were constructed using 
the older OPAL equation of state, which did not treat relativistic effects at temperatures 
and densities relevant to the solar core. This resulted in a somewhat deficient core structure 
(see EUiot & Kosovichev, 1998). We shall discuss in more detail below results obtained from 
models that were instead constructed with the corrected OPAL equation of state. 

Model S06-(-(AGS, OP) was made because a localized change in opacities can resolve 
the problem of the convection zone depth in models with AGS abundances (see, e.g. Basu 
& Antia 2004; Bahcall et al. 2005a). However, it is very clear from the fine spacings, and 
separation ratios that this model fails in the core. 

We find that the model BBS05-3 does reasonably. This model was constructed by 
increasing the amount of neon compared to the AGS abundances, in addition to increasing 
the quantities of some of the other elements. A strategy of increasing the amounts of neon 
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and other elements (within their error bars) was suggested by Antia & Basu (2005) and 
Bahcall et al. (2005) as a potential way to solve the disagreement between the helioseismic 
observations and the models with AGS abundances. Neon contributes to the opacity near the 
convection-zone base. An increase in the neon abundance therefore deepens the convection 
zone, bringing the models into agreement with the helioseismically determined convection- 
zone depth. It is worth adding that because of the abundance increases in some of the 
elements this model has a relatively high ZjX (0.02069). 

Differences between the results of the various models ultimately have to be understood 
in terms of the sound-speed gradient in the core which, as per Eq. 2, determines the fine 
spacings and therefore also the separation ratios. Sound speed is determined both by tem- 
perature T and mean molecular weight /x, with oc Tf/i. The gradient of sound speed is 
therefore determined by the gradients of temperature and /i as well. Although the sound- 
speed gradient is negative in the bulk of the Sun, it is positive at very small radii because of 
the different radial dependences of T and with the increase in /i at smaller r compensating 
the increasing T to give a lower c. 

We can use the three YREC models to determine why the other models with lower 
metallicity fare worse than models with higher metallicity in comparison with the BiSON 
data. Fig. 8 shows the fine spacings and separation ratios of the YREC models, together 
with the BiSON data (points with error bars). We can see that increasing the surface Z/X 
decreases the fine spacings and separation ratios. To explain these results we show in Fig. 9 
the sound-speed gradient and the mean molecular weight profiles of the YREC models. 
It is clear that differences between the models lie very close to the center, where dc/dr 
is positive. A more detailed analysis of the structure of the YREC models shows that the 
term in the sound-speed gradient is instrumental in causing the differences, 
in particular the factor y^Yfjj^. The models have very similar temperatures, temperature 
gradients and /i gradients, but the larger /i in models with higher Z/X increases dc/dr in 
the core, thereby decreasing the value of the integral in Eq. 2. 

The jJL dependence of the sound-speed gradient explains why models BP04-I- and BBS05- 
3 do not behave like BS05(AGS, OP) and S06+(AGS,OP). Model BP04+ does not have as 
low metallicities as BS05(AGS, OP) and S06-|-(AGS,OP) and hence fares better than these 
models. That said, the metallicity of BP04-I- is low enough to give larger fine spacings and 
separation ratios than are seen in the GS models. 

Model S06-|-(AGS, OP) does not do well in the core, despite satisfying the helioseismic 
constraints at the base of the convection zone, because /x in its core remains similar to that 
in model BS05(AGS, OP). The change in opacity in S06+(AGS, OP) does however induce 
changes in the rest of the model (solar structure is, after all, determined by a set of highly 



- 12 - 



coupled equations), and this causes enough difference in the core temperature gradient to 
give different fine spacings and separation ratios compared with BS05(AGS, OP). 

This leads us to the question of the differences in the spacings and ratios of models 
BP04 and BS05(OP). That there are differences is not surprising because even though they 
have the same metaUicity, the models were constructed with shghtly different opacities. The 

temperature gradient in the radiative zone (and that includes the solar core) is determined 
by the opacity, and hence different opacities will imply a different temperature gradient. 
Differences in the treatment of diffusion and ^^O burning also have a small effect because 
both change the /i gradient slightly. Nonetheless, the overall opacity differences in the models 
lead to only modest differences in the fine spacings and separation ratios. 

From the above discussion it is clear that models with low metallicities, such as those 
with AGS abundances, do not satisfy the observational constraints imposed by the BiSON 
fine-spacing and separation-ratio data. Though subtle changes in physics (such as in the 
formulation of diffusion, etc.) do change the core structure, the changes are sufficiently 
small to allow us to distinguish between the GS and AGS abundances. In fact, looking at 
Fig. 8, we might be tempted to determine the solar metaUicity using the fine spacings and 
separation ratios. 

4. Implications for neutrino uncertainties 

Uncertainties in solar neutrino predictions arise from uncertainties in inputs to the solar 
models. It has, however, been long recognized that the predominant source of uncertainty 
is the composition (Sears 1964, Bahcall 1966). As a result there has been a large body of 
work devoted to evaluating neutrino uncertainties, among the more recent being Fiorentini & 
Ricci (2002), Couvidat ct al. (2003), Boothroyd & Sackmann (2003), Bahcall & Pinsonneault 
(2004), Young & Arnett (2005), Bahcall & SereneUi (2005), and Bahcall et al. (2006). 

Bahcall & SereneUi (2005) performed a very detailed and rigorous investigation of how 
the abundances of individual elements affect solar neutrino predictions. Their study was 
motivated by the fact that models with AGS abundances show discrepancies with respect 
to the Sun just below the solar convection zone. Bahcall & SereneUi (2005) found that the 
largest uncertainties in the ^Be and neutrinos are due to the uncertainty of the solar 
iron abundance. However, uncertainties in oxygen, neon, silicon and sulphur also contribute 
significantly to the ^Be and ^B fiux uncertainties. The p-p neutrino fiux is most sensitive the 
changes in the iron abundance, but because carbon uncertainties are larger, both elements 
dominate and contribute comparable amounts to the total uncertainty. It should be noted 



-13- 



that the neutrinos are the ones detected by the water (Kamiokande & Superkamiokande) 
and heavy-water detectors (Sudbury Neutrino Observatory). ^B, ^ Be and pep neutrinos are 
detected by Chlorine detectors. In order to calculate the uncertainty, Bahcall & Serenelli 
(2005) adopted what they called a "conservative" uncertainty in the abundances. This 
uncertainty is basically the difference between the GS and AGS abundance values. Their so- 
called "optimistic" uncertainty was based on the error bars on the abundance from the AGS 
table. The "conservative" uncertainties are much larger than the "optimistic" uncertainties. 
It should be noted that quoted uncertainties in the GS and AGS tables are fairly similar. 
Work done by Bahcall et al. (2006) confirmed these results through an elaborate Monte-Carlo 
simulation. 

Our current work has shown that models with AGS abundances fare much worse than 
models with GS abundances. We also see, from Fig. 8, that the uncertainties in the abun- 
dances are probably much less than the difference between the GS and AGS values. We 
therefore conclude that the "conservative" uncertainties in the neutrino fluxes are unduly 
pessimistic and can be lowered substantially. Not having done the full analysis, we are un- 
able to say if we are justified in lowering the levels to the "optimistic" levels of Bahcall & 
Serenelh (2005) and Bahcall et al. (2006). 

Our conclusion however comes with one condition, which is related to the model with 
modified neon abundance (BBS05-3) and the fact it did well in comparisons with the BiSON 
data. The solar neon abundance is very uncertain since the abundance cannot be measured 
at the photosphere. This had led Antia & Basu (2005) and Bahcall et al. (2005b) to propose 
increasing the neon abundance as a solution to the helioseismic problems posed by the AGS 
models. Drake & Testa (2005) found that most neighboring stars seem to have a much 
higher Ne/0 ratio compared to the Sun, which supported increasing the neon abundance. 
However, Schmelz et al. (2005) and Young (2005) each reanalyzed solar X-ray and UV data 
and found that the Ne/0 ratio of the Sun is indeed low. Bochsler et al. (2006) also suggest 
a low solar Ne/0 abundance. A more detailed analysis of the convection zone also points to 
the fact that an enhanced neon abundance will not reconcile models with the helioseismic 
data (Basu & Antia 2006). Thus although the neon issue has not been resolved completely, 
it looks increasingly unlikely that the solar neon to oxygen ratio is high. We are therefore 
reasonably confident about our claim that uncertainties in the predicted solar neutrino fiuxes 
can be lowered. 
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5. 



Conclusions 



We have used BiSON data collected over 4752 days to determine the frequencies of low- 
degree acoustic solar oscillation modes. These frequencies have been used to calculate the 
fine or small frequency spacings and frequency separation ratios for the Sun. These spacings 
and ratios have then been used to test different solar models to investigate whether or not 
we can constrain solar abundances. 

We find that models constructed with the older Grevesse & Sauval (1998) mixture satisfy 
the fine-spacing and separation-ratio constraints much better than do models with the newer 
Asplund et al. (2005) abundances. In fact, models with high metaUicity constructed with 
outdated opacities and equation-of-state data also fare better than do the recent models 
with Asplund et al. (2005) abundances. Our investigation shows that the fine spacings and 
separation ratios depend sensitively on the metaUicity, and low-metallicity models can be 
ruled out. 

Our results lead us to conclude that uncertainties on predicted solar neutrinos, which 
had been raised because of the Asplund et al. (2005) abundances, can now be lowered. 

This paper utihzes data collected by the Birmingham Solar- Oscillations Network (Bi- 
SON), which is funded by the UK Particle Physics and Astronomy Research Council (PPARC) 
We thank the members of the BiSON team, colleagues at our host institutes, and all others, 
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of PPARC. The radio flux observations are made at Penticton by the National Research 
Council of Canada and are available from the World Data Center. SB acknowledges partial 
support from NSF grant ATM-0348837. She would also like to thank the BiSON group for 
their hospitality during the time this work was carried out. AMS is partially supported by 
the NSF (grant PHY-0503684), the Association of Members of the Institute for Advanced 
Study and the W. M. Keck Foundation through a grant-in-aid to the Institute for Advanced 
Study. 
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Table 1. Difference in fine-structure spacings (BiSON minus models) 







< J 02(11 


) 
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Model 


Weighted mean 


Weighted 


Weighted mean 


Weighted 




difference (/xHz) 


RMS (/xHz) 


difference (/xHz) 


RMS (/xHz) 


JCDS 


-0.058 


(-10.8(7) 


0.072 (13.3(7) 


-0.033 


(-4.5a) 


0.053 (7.2a) 


SAC 


-0.133 


(-24.7(j) 


0.139 (25.9(7) 


-0.083 


(-11.2a) 


0.094 (12.8a) 


BP04 


-0.009 


(-1.7(7) 


0.047 (8.7(7) 


-0.043 


(-5.8a) 


0.091 (12.3a) 


BP04+ 


-0.102 


(-18.9(7) 


0.126 (23.4(7) 


-0.204 


(-27.7a) 


0.237 (32.2a) 


BS05(OP) 


-0.051 


{-9.5a) 


0.065 (12.0(7) 


-0.083 


(-11.3a) 


0.114 (15.5a) 


BS05(AGS,OP) 


-0.230 


(-42.6(7) 


0.243 (45.1(7) 


-0.360 


(-48.7a) 


0.382 (51.7a) 


BBS05-3 


-0.034 


(-6.2a) 


0.063 (11.7a) 


-0.078 


(-10.5a) 


0.136 (18.4a) 


S06+(AGS,OP) 


-0.145 


(-27.0(7) 


0.160 (29.7a) 


-0.200 


(-27.1a) 


0.267 (36.2a) 



Table 2. Difference in frequency separation ratios (BiSON minus models) 



ro2{n) nsin) 



Model Weighted mean Weighted Weighted mean Weighted 

difference (%) RMS (%) difference (%) RMS (%) 

JCDS -0.43 (-9.5a) 0.63 (13.7a) -0.03 (-0.8a) 0.23 (5.9a) 

SAC -0.99 (-21.6a) 1.08 (23.5a) -0.22 (-5.6a) 0.29 (7.5a) 

BP04 +0.03 (+0.6a) 0.26 (5.7a) +0.02 (+0.4a) 0.17 (4.5a) 

BP04+ -0.85 (-18.6a) 0.92 (20.0a) -0.84 (-22.1a) 0.89 (23.2a) 

BS05(OP) -0.39 (-8.5a) 0.46 (10.0a) -0.29 (-7.5a) 0.35 (9.2a) 

BS05(AGS,OP) -2.02 (-43.9a) 2.06 (44.8a) -1.76 (-45.8a) 1.80 (46.9a) 

BBS05-3 -0.05 (-1.1a) 0.27 (5.9a) +0.02 (+0.4a) 0.22 (5.7a) 

S06+(AGS,OP) -0.98 (-21.3a) 1.04 (22.5a) -0.58 (-15.1a) 0.79 (20.6a) 
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Fig. 1. — Differences between fine spacings made with and without the solar-cycle correction 
(in the sense corrected spacings minus raw spacings). Left-hand panel: difference/correction 
in do2{n). Right-hand panel: difference/correction in ^13(71). 
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Fig. 2. — Solid line in each panel: fine spacings do2{n) of solar model identified in plot title. 
Points with error bars: (io2(^) from the corrected BiSON frequencies. 
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Fig. 3. — Solid line in each panel: fine spacings (ii3(n) of solar model identified in plot title. 
Points with error bars: diz{n) from the corrected BiSON frequencies. 
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Fig. 4. — Solid line in each panel: separation ratios rQ2{'n) of solar model identified in 
title. Points with error bars: ro2(n) from the corrected BiSON frequencies. 
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Fig. 5. — Solid line in each panel: separation ratios ri3(n) of solar model identified in 
title. Points with error bars: r 13(71) from the corrected BiSON frequencies. 
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Fig. 6. — Residuals given by subtraction of best-fit straight line (to fine spacings ver- 
sus frequency) from BiSON fine spacings and each set of model fine spacings. Left-hand 
panel: residuals in do2{n). Right-hand panel: residuals in di3{n). BiSON data are plotted 
with their associated error bars. Various models rendered as follows: JCDS (red, solid); 
SAC (red, dashed); BP04 (cyan, sohd); BP04+ (black, sohd); BS05(OP) (cyan, dashed); 
BS05(AGS,OP) (black, dotted); BBS05-3 (green); S06+(AGS,OP) (black, dashed). 




Fig. 7. — Differences between observed BiSON and model separation ratios (in sense BiSON 
minus model). Left-hand panel: differences in ro2(?T-). Right-hand panel: differences in ri3(r2). 
Various models rendered as follows: JCDS (red, solid); SAC (red, dashed); BP04 (cyan, 
solid); BP04+ (black, solid); BS05(OP) (cyan, dashed); BS05(AGS,OP) (black, dotted); 
BBS05-3 (green); S06+(AGS,OP) (black, dashed). 
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Fig. 8. — Impact of changes to Z/ X in the models. The upper panels show the fine spacings of 
three models computed with the YREC evolution code and the OPAL equation of state, but 
with different Z/X. Solid hne: model YREC(GS), with Z/X = 0.0229. Dotted fine: model 
YREC(AGS), with Z/X = 0.0165. Dashed line: model YREC(0.03), with Z/X = 0.03. 
Lower panels: same, but for the separation ratios. 
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Fig. 9. — Left-hand panel: The function {l/r){dc/dr) for the three YREC models. This 
function determines to a large extent the fine spacings, and therefore also the separation 
ratios. Right-hand panel: The mean-molecular weight profile of the three YREC models. 



